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Abstract. 

In recent years, X-ray observations with Chandra and XMM-Newton iiave significantly increased our understanding of 
rotation-powered (radio) millisecond pulsars (MSPs). Deep Chandra studies of several globular clusters have detected X-ray 
counterparts to a host of MSPs, including 19 in 47 Tuc alone. These surveys have revealed that most MSPs exhibit thermal 
emission from their heated magnetic polar caps. Realistic models of this thermal X-ray emission have provided important 
insight into the basic physics of pulsars and neutron stars. In addition, intrabinary shock X-ray radiation observed in "black- 
widow" and peculiar globular cluster "exchanged" binary MSPs give interesting insight into MSP winds and relativistic shock. 
Thus, the X-ray band contains valuable information regarding the basic properties of MSPs that are not accesible by radio 
timing observations. 
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INTRODUCTION 

Rotation-powered ("recycled") millisecond pulsars 
(MSPs) were discovered at X-ray energies during the 
ROSAT all-sky survey iH. In its lifetime, ROSAT de- 
tected X-ray emission from 10 MSPs in the field of 
the Galaxy and one in the globular cluster M28 ij]. 
These observations established that MSPs are very faint 
X-ray sources with typical X-ray luminosities in the 
range Lx ^ 10^"^^^ ergs s^'. In recent years, the superb 
sensitivity of the Chandr^ and XMM-Newtof^ X-ray 
observatories have permitted much more detailed studies 
of these objects. Chandra, in particular, with its unprece- 
dented sub-arcsecond angular resolution, has proven to 
be ideally suited for studies of globular clusters, where 
MSPs are found in great abundance. In addition, both 
facilities have allowed detailed spectroscopic and timing 
studies of the nearby MSPs in the field of the Galaxy 
Bj, 3] . Herein, we summarize the current state of affairs 
and examine future prospects in the X-ray study of 
MSPs. 



X-RAY STUDIES OF GLOBULAR 
CLUSTER MILLISECOND PULSARS 

Radio timing surveys have revealed that globular clus- 
ters are veritable treasure troves of MSPs [e.g.,|5llal- As 
such, they allow a very efficient way to study MSPs at 
X-ray energies due to their inherent faintness of these 



objects, which necessitates very long exposures. This is 
best examplified by the detection of X-ray counterparts 
to all 19 MSPs with known radio timing positions (Fig. 
1) in a deep Chandra ACIS-S exposure of 47 Tuc. The 
large majority of the 47 Tuc MSPs (16 out of 19) ap- 
pear to be soft, thermal X-ray sources with temperatures 
TeS ^ (1 — 3) X 10^ K, emission radii Reff^ 0.1—3 km, 
and luminosities Lx '^ 10^"^^^ ergs s^' JTll. The inferred 
effective emission areas imply that this radiation orig- 
inates from the heated magnetic polar caps of the pul- 
sar, as suggested by theoretical pulsar models |8, 9]. The 
thermal nature of the observed X-ray emission was sub- 
sequently confirmed with Chandra HRC-S timing obser- 
vations, which revealed relatively broad pulsations with 
relatively low (<50%) pulsed fractions 1.10.1 . For sev- 
eral MSPs, there is indication of an additional thermal 
component (implying non-uniform heating of the po- 
lar cap), similar to what is seen in several nearby field 
MSPs yl]. These similarities in X-ray properties imply 
that, in general, no systematic differences exist between 
MSPs in globular clusters and those in the field of the 
Galaxy. It has been speculated that globular cluster MSPs 
should differ from those in the field of the Galaxy due 
to repeated dynamical binary exchanges and accretion 
episodes, which may alter the magnetic field topo log y 
of the pulsar from a dipolar to a multipolar one Illll[l2ll . 
However, the X-ray data indicates that even if most glob- 
ular cluster MSPs do undergo "re-recycling", this does 
not significantly alter the observable properties of these 
objects. 
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FIGURE 1. Chandra ACIS-S 0.3-6 keV image of the core of the globular cluster 47 Tucanae. The 1" circles are centered on the 
radio timing positions of the 19 MSPs. An X-ray source is associated with the position of each MSP [see 7, for details]. 



MODELING THERMAL X-RAY 

EMISSION FROM MILLISECOND 

PULSARS 

It has been shown that modeUng X-ray spectroscopic and 
timing data of thermal MSPs can be used to extract im- 
portant information regarding the properties of the neu- 
tron star [IJ, [3. uM- Given that the vast majority of 
MSPs appear to be thermal X-ray sources, it is impor- 
tant to develop realistic theoretical models to ascertain 
whether this method can be employed to gain useful in- 
sight into the properties of MSPs. Such models of surface 
emission from MSPs need to incorporate general rela- 
tivistic effects of gravitational redshift and bending of 
photon trajectories as well as special relativistic effects 
(Doppler boosting and abeiTation) due to the rapid mo- 
tion of the stellar surface. It is also necessary to consider 
an optic ally- thick, weakly magnetized light-element at- 
mosphere [see, e.g., [T6|]. Such an atmospheric layer is 
expected at the surface of a MSP given the standard for- 
mation scenaro of MSPs, involving accretion of a sub- 
stantial amount of gas [1711. 

The synthetic pulse profiles of MSPs generated using 
these models are characterized by broad pulsations and 
moderate pulsed fractions, which are significantly larger 
than those of a blackbody emission spectrum for the 
same assumed parameters (see Fig. 2). Most notably, 
unlike the beamed radio emission, the duty cycle of the 
X-ray radiation is always 100% and is observable for 



all combinations of a and ^. These characteristics have 
important implications for population studies of MSPs 
(see below). 

A hydrogen atmosphere polar cap X-ray emission 
model has been found to be in fair agreement with 
the available X-ray data of the three nearest known 
MSPs, PSRs 10437^715 H [H, J2124-3358, and 
J0030H-0451 1 19]. On the other hand, a blackbody emis- 
sion model is inconsistent the observed pulse profiles. 
This provides compelling evidence that a gaseous atmo- 
sphere is indeed present at the surface of these neutron 
stars. For PSRs J0437^715 and J2124-3358, the data 
indicates that the two polar caps are not diametrically 
opposite on the stellar surface, implying an off-center 
dipole magnetic axis. Perhaps more importantly, the ther- 
mal X-ray pulsations offers interesting constraints on the 
mass-to-radius ratio {M/R) of the neutron star In partic- 
ular, for J0437-4715, J2124-3358, and J0030H-0451 we 
find R = 6.9- 10.6, R > 7.2, and R > 9.5 km, respec- 
tively (68% confidence), assuming M — 1.4 M0. 

Modeling of the thermal X-ray pulse profiles from 
MSPs also gives useful limits on the magnetic field ge- 
ometry of the pulsar, which in turn provides stringest 
observational constraints on theoretical pulsar models. 
For example, Ruderman has postulated the presence of 
crustal "plates" on the NS surface formed by shear 
stresses on the crust caused by neutron superfluid vor- 
tex lines pinned to lattice nuclei. The motion of these 
plates would cause the magnetic fields of MSPs to mi- 
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FIGURE 2. Representative hydrogen atmosphere model 
lightcurves for a rotating M = 1.4 Mq, i? = 10 km NS with 
two point-like antipodal hot spots [18]. The panels show pulse 
profiles for (a = 10°, f = 30°), (a = 30°, ^ = 60°), (a = 60°, 
^ = 80°), and (a = 20°, ^ = 80°), from top to bottom, respec- 
tively. The solid lines in each plot correspond to a H atmosphere 
(black) and blackbody emission (green) with no Doppler effect 
included. The dashed lines show the effect of Doppler boosting 
and aberration for P = 4 ms. All fluxes are normalized to the 
value for a = ^ = 0. 



grate across the stellar surface, resulting in either an 
aligned magnetic field or one pinched at one of the spin 
poles. However, the observed thermal pulse shapes and 
pulsed fractions from the nearest MSPs (~30-50%) sug- 
gest that the magnetic field resembles a (nearly) cen- 
tered dipole and is moderately inclined with respect to 
the spin axis (see Fig. 3). This, in turn, implies that the 
dipole fields of these (and likely all) MSPs do not have 
a tendency to migrate towards the spin pole or align 
with the spin axis. Moreover, the infered X-ray-emitting 
areas are consistent with the classical polar cap area 
[Rp^ = {InR/cPyl^R km] and do not favor so-called 
"squeezed" polar caps ll20ll . 



"BLACK WIDOW" AND "EXCHANGED" 
MSP BINARIES 

A number of radio eclipsing binary MSPs, such as the 
canonical "black widow" pulsar B1957H-20 iMIll, 47 



Tuc J and O 101, and M71A m, exhibit unresolved hard 
non-thermal X-ray emission. This is likely synchrotron 
radiation from a shock formed by the action of the pulsar 
wind on the companion star [24]. X-ray studies (comple- 
mented by optical observations) of these systems offer a 
useful diagnostic of the properties of MSP winds. 

The eclipsing binary MSPs PSR J1740-5340 in the 
globular cluster NGC 6397 1251 [lllj and PSR J0024- 
7204W in 47 Tuc 1.26.12711 exhibit very peculiar radio, op- 
tical, and X-ray properties that set them apart from other 
MSPs (see Fig. 4). The unusual behavior of these systems 
most probably stems from the fact that they are bound 
to unevolved (main-sequence-Uke) companions. This is 
at odds with the standard formation scenario of MSPs, 
which requires the companion to evolve off the main se- 
quence in order for accretion and spin-up to commence 
lUTIl . It is likely that the current binary partners were not 
the stars that spun up the MSPs. As the dense environ- 
ment of a globular cluster core is conducive to close dy- 
namical encounters between stars, it is very likely that 
the present companion was acquired in a close binary- 
binary exchange encounter. These "exchanged" binaries 
offer a unique opportunity to study NSs transitioning 
from accretion to rotation power. In particular, the ob- 
served X-ray and optical properties of the radio MSP 
PSR J0024-7204W are remarkably similar to those of 
the low-mass X-ray binary and X-ray millisecond pul- 
sar SAX J1808. 4-3658 in quiescence. This supports the 
conjecture that the non-thermal X-ray emission and op- 
tical modulations seen in the SAX J1808.4-3658 sys- 
tem in a quiescent state are due to interaction between 
the wind from a reactivated rotation-powered pulsar and 
matter from the companion star. The striking similarities 
between the two systems provide further support for the 
long-sought connection between millisecond radio pul- 
sars and accreting neutron star systems. 



NON-THERMAL X-RAY PULSATIONS 

A very small subset of MSPs show remarkably sharp X- 
ray pulsations, indicative of non-thermal magnetospheric 
emission processes. These MSP include the exception- 
ally energetic PSRs B1937H-21 0], B1821-24 [29], and 
J0218H-4232 {E = 2.4 x 10^5 ergs s"', B^„,-f = 4.3 x 10** 
G) [28] as well as the mildly recycled relativistic binary 
MSP PSR J0737-3039A |30]. This striking difference in 
X-ray spectral properties compared to most other MSPs 
implies a substantially different physical regime in the 
pulsar magnetosphere. Models of pulsars suggest that the 
combination of spin properties of these four pulsars place 
them near the locus of theoretical curvature radiation 
(CR) pair-production death lines [31] on the P — P di- 
agram. Therefore, in these of pulsars, pair-production in 
the magnetosphere probably occurs via curvature radia- 
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FIGURE 3. {Left) XMM-Newton EPIC-pn 0.3-2 keV pulse profile of PSR J0437-4715. The blue line shows the best fit curve for 
an off-center dipole field and a H atmosphere, while the green line is of a blackbody that best fits the spectrum of J0437-4715 for 
the same assumed compactness and geometry. The red line corresponds to the best fit H atmosphere model for a centered dipole 
(Right) Equal-area sinusoidal projection surface emission map of PSR J0437-4715 as inferred from the X-ray data for R = 10 km 
and M = 1.4 Mq. The dashed line shows the expected antipodal position of the secondary polar cap. 



in 0.002 



> 
<u 

>, 

QD 
1- 
<U 

c 

c 
o 

o 

s: 



0.001 



6000 



4000 



2000 ^,.-j 



—\ — I — I — r- 




^••^V.. •••=•• .V ■V""^v.. •••••■ .V 

I I I I I I I I I I I I I I I I L_ 



%'si^ 



0.5 1 1.5 2 

Orbital phase 

FIGURE 4. X-ray lightcurve of 47 Tuc W folded at the 
binary period (Pj, = 3.2 h). Note the remarkable absence of 
hard X-ray photons between phases 0.4 and 0.7, which can 
be explained by a geometric occultation of the intrabinary 
shock by the unusual main-sequence secondary star [see |23, 
for details]. 



tion emitted by accelerating electrons. On the other hand, 
most MSPs (such as those in 47 Tuc), lie well below the 
theoretical CR pair production death lines. In these ob- 
jects, inverse Compton scattering (ICS) is the primary e^ 
production mechanism and the conditions in the magne- 



tosphere favor substantial heating of the polar caps by a 
return flow of e+ istl . 



FUTURE PROSPECTS 

Due to the inherent faintness of MSPs, at present, 
only a handful of nearby objects can be studied in 
great detail [3j. The next generation of X-ray facilities. 
Const ellation-X and XEUS, are expected to have ~100 
times greater sensitivity than which would enable stud- 
ies of a much greater sample of MSPs. More impor- 
tantly, these telescopes could be employed in blind tim- 
ing searches for nearby (<2 kpc) MSPs that are not ob- 
servable at radio frequencies due to unfavorable view- 
ing geometry. This is possible because the effect of light 
bending combined with the (nearly) antipodal configu- 
ration of the two MSP hot spots ensure that the thermal 
radiation is seen by all distant observers for any com- 
bination of a and ^. In contrast, at radio frequencies a 
pulsar is not observable if |a — ^| exceeds the opening 
half-angle p of the radio emission cone. 

Figure 5 shows the X-ray pulsed fraction of a 10 km, 
1.4 Mpj MSP as a function of a and ^. Also shown are 
lines delineating the region of the a — ^ plane for which a 
pulsar with a given radio emission cone width is observ- 
able at radio frequencies. If we assume a uniform distri- 
bution of pulsar obliquities (a) and viewing angles (Q, 
for p ~ 30° a substantial portion (--45%) of the MSP 
population is invisible to us in the radio. On the other 
hand, if we consider an X-ray timing survey with a limit- 
ing pulsed fraction sensitivity of ^10%, only ^^5-15% 
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FIGURE 5. The a - ij plane for MSPs. The solid lines 
represent contours of constant thermal X-ray pulsed fraction 
for i? = 10 km and M = 1.4 Mj. The dot-dashed, dashed, and 
dotted lines show 10°, 30°, and 60°, respectively, opening half- 
angles of the radio emission cone. 



(depending on the mass and radius of the star) of the 
MSPs will go undetected as pulsed sources though they 
will still be detected as X-ray sources fl9j. The Galactic 
population of MSPs may in fact be preferentially clus- 
tered in a certain range of a due to the poorly understood 
effects of the accretion and magnetic field reduction pro- 
cesses during the LMXB phase on the NS. A deep X-ray 
timing survey of nearby (<l-2 kpc) MSPs may, in prin- 
ciple, reveal whether this is indeed the case. 



CONCLUSIONS 

To date, nearly 40 MSPs have been detected at X-ray 
energies, 19 of which are found in the globular cluster 
47 Tuc \^. It has become apparent that X-ray spectro- 
scopic and timing observations of MSPs represent a very 
promising approach towards elucidating important pul- 
sar properties that are not measurable by other observa- 
tional means such as radio timing observations. Further 
study of these sources will provide more important in- 
formation regarding the X-ray properties of the Galactic 
population of MSPs, as well as valuable insight into the 
fundamental properties of NSs. 
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